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Variational bounds for the shear viscosity of gelling melts 
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PACS 64 . 60 . Ht - Dynamic critical phenomena 

PACS 61.25.Hq - Macromolecular and polymer solutions; polymer melts; swelling 
PACS 61.20.Lc - Time-dependent properties; relaxation 

Abstract. - We study shear stress relaxation for a gelling melt of randomly crosslinked, interacting 
monomers. We derive a lower bound for the static shear viscosity n, which implies that it diverges 
algebraically with a critical exponent k 2v — f3. Here, v and (5 are the critical exponents of 
percolation theory for the correlation length and the gel fraction. In particular, the divergence is 
stronger than in the Rouse model, proving the relevance of excluded-volume interactions for the 
dynamic critical behaviour at the gel transition. Precisely at the critical point, our exact results 
imply a Mark-Houwink relation for the shear viscosity of isolated clusters of fixed size. 
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Introduction. — The most prominent feature of a 
gelling (macro-) molecular liquid is the divergence of the 
static shear viscosity as the gelation transition is ap- 
proached from the sol side. The melt becomes increas- 
ingly more viscous, reflecting larger and larger solid like 
regions and signalling the emergence of shear rigidity in 
the amorphous solid state of the gel. Many experiments 
have tried to measure the divergence of the shear viscosity 
quantitatively, yet the data are still very confusing. 

The gelation transition is known to be a percolation 
transition, which is controlled by the crosslink density c, 
see e.g. [1] for a review. It has been analysed with the 
tools of critical phenomena, in particular it has been ver- 
ified that the cluster size distribution follows the scaling 
law f(n) ~ n~ T exp{— n/n*} if the distance to the critical 
point, c cl it, is sufficiently small. Here n is the cluster size 
or more precisely the number of crosslinked monomers of a 
cluster. The quantity n* ~ |c— c cr i t |~ 1 / CT has the meaning 
of a "cut-off size" and a and r are critical exponents [2] . 

The divergence of the static shear viscosity is usu- 
ally analysed in terms of an algebraic singularity r\ ~ 
(ccrit — c)~ k and a dynamic critical exponent k. However 
the values of k which have been deduced from experiment 
vary enormously - in contradiction to the concept of uni- 
versal exponents put forward by the dynamic renormalisa- 
tion group theory. To appreciate the scatter of the data we 
summarise the reported values for k in Table 1. We have 
also included results from simulations in Table 1, based 
either on the bond-fluctuation model [3] or on molecular 



dynamics in combination with a soft pair potential [4,5]. 

The origin of the broad scatter of the data is not under- 
stood at present. It has been suggested that there might 
be more than one universality class for the dynamic critical 
exponents. For example, in dense solutions hydrodynamic 
interactions play a crucial role, whereas in melts such in- 
teractions are absent. However, it is not obvious, how to 
group the above data into two universality classes. An- 
other explanation for the scatter has evoked the size of 
the critical region, which shrinks with growing linear size 
of the crosslinked (macro-) molecules [19]. Consequently 
some experiments could be outside the critical region or 
measure crossover exponents. 

In view of the controversial results it is highly desirable 
to obtain rigorous results for the divergence of the static 
shear viscosity. In this paper we consider a Rouse model, 
generalised to account for excluded-volume interactions, 
and derive an exact lower bound for rj. This allows us to 
conclude that the viscosity diverges at least as strongly 
as the lower bound, implying the inequality k 2v — (3 
for the critical exponent k. Here, v and (3 are the crit- 
ical exponents of percolation theory for the correlation 
length and the gel fraction. In d = 3 spatial dimensions 
the lower bound for the exponent implies k 1.32, in dis- 
agreement with several entries in Table 1. Furthermore the 
lower bound exceeds the value for the critical exponent of 
the Rouse model without excluded volume, proving that 
excluded-volume interactions are indeed relevant for the 
dynamic critical behaviour. We also derive a lower bound 



p-1 



C. H. Kohler, H. Lowe, P. Miiller, A. Zippelius 



Table 1: Numerical values for the viscosity exponent from experiments and simulations 
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for the case of a single (isolated) cluster with n monomers. 
This lower bound is found to agree with an upper bound, 
resulting in an exact expression for the shear viscosity rj n ^ s 
for single clusters of fixed size n. If the radius of gyration 
scales with cluster size according to i? gyr oc n 1 / dl - e , where 
d{. s is the Hausdorff fractal dimension for isolated clus- 
ters, then the single-cluster viscosity is shown to obey the 
Mark-Houwink scaling relation rj n ^ s ~ n 2 / df - s . 

Model. — We consider a system of N identical par- 
ticles at positions {R\, . . . ,Rn} in a <i-dimensional vol- 
ume V. Permanent crosslinks connect M randomly chosen 
pairs of particles J\f = {(h,ji), {12,32), -(im, 3m)}- The 
crosslinks are modelled by harmonic springs 



1 M 

e=l 



(i) 



of typical length a. (Energies are measured in units of 
k^T). A particular crosslink realisation J\f can be decom- 
posed into L clusters of crosslinked particles: M = \jf =1 J\fi . 
The Z-th cluster, J\fi, contains Ni particles. All the parti- 
cles interact via an isotropic pair potential 



C/ EV := u (\ R i- R A) 



(2) 



which can be taken as a Lennard- Jones potential or just 
a as hard core to ensure excluded volume. The details 
of the potential are irrelevant for this paper. The total 
potential energy is the sum of the harmonic crosslinks and 
the excluded volume: U := U x + Uev ■ 

Shear viscosity. — We are interested in shear flow 
and impose an external velocity field which points in the 
^-direction and increases linearly with y: v(r) := nye x . 
Here k denotes the shear rate which is taken time inde- 
pendent because we aim to compute the static shear vis- 
cosity. In the following we shall need the symmetrised ex- 
ternal flow field < ym := kK^v^ with K vtt := (S^S^y + 
&v,y&n,x) I?* and summation over repeated Greek indices 
implied. The monomers of the melt perform a Brownian 
motion on a semi-microscopic scale. The monomer config- 
uration at time t is then described by the A-particle prob- 
ability density ip(Ri, i?2, Rn, t), which obeys the time- 
dependent Fokkcr-Planck or Smoluchowski equation [20] 

Here £ > denotes the friction constant. 



The external velocity field builds up stresses in the sys- 
tem. Following Kirkwood, they are described by the stress 
tensor er with components [20] 



VjJ, 



1 



dU 



(4) 



The stress tensor is proportional to the monomer density 
p, and the brackets (-} stat indicate the thermal average 
over particle positions with respect to the suitably nor- 
malised stationary solution V'stat of (3) in the presence of 
a shear flow. Note that the stress tensor as defined in (4) 
is symmetric. The intrinsic linear-response viscosity is ob- 
tained from the stress tensor (4) as 77 := lim K ^ & xy / 

In order to determine the shear viscosity rj within linear 
response, we follow [21] and expand the stationary proba- 
bility density V'stat = V'cq(l + K<fi) + 0(k 2 ) in powers of the 
shear rate. The zeroth-order term is given by the canonical 
equilibrium state V'cq ~ exp{— U}. In order to ensure the 
correct normalisation of V'stat we require (</>) = 0, where 
(•) denotes the thermal average over {Ri} weighted by 
V'eq- The above expansion allows us to express the linear 
viscosity as an average over the equilibrium distribution 
in two alternative forms 



OR? 7 (N^ i \dR l <)R, /, 



eq 



oq 

(5) 

Here only the symmetrised external flow field enters, even 
if the externally imposed flow is not symmetric. This is a 
direct consequence of the symmetry of the stress tensor as 
defined in (4). 

All the difficulties in the computation of n are hidden 
in the stationary non-equilibrium distribution <f>. In this 
paper we shall follow Fixman [21] and derive variational 
bounds for r\ which circumvent the exact knowledge of 
4>. Instead a guess for <ft will serve as a test function in 
the variational principle with the quality of the bounds 
depending strongly on the chosen test function. To get a 
first hint, we consider the simpler case without excluded- 
volume interactions ("phantom clusters"), which can be 
solved exactly for the viscosity. In this case the centre of 
mass positions of the clusters 



Ni Z - 



(6) 



follow the external flow field, giving rise to a nonzero 
current in the stationary state. The stationary distribu- 
tion Vstat depends only on the particle positions relative 
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to the centre of mass of the cluster to which they be- 
long, i.e. on SRi := Ri — R^ l \ It is explicitly given by 
V'stat ~ (1 + K<p x ) exp{-J7 x } with 



1=1 ieM 



(7) 



Inserting (7) into (5) and using isotropy, we recover the 
result 



i=i 



(8) 



for the viscosity of the Rouse model for gelling polymers 
[22-24] in terms of the cluster-weighted sum of the squared 
radii of gyration 



(9) 



Bounds for the viscosity. — To derive a lower bound 
for the viscosity in the presence of excluded- volume inter- 
actions, we define the functional 

(10) 

for smooth scalar fields 4>. With help of the two alter- 
native expressions for the viscosity (5), the latter can be 
represented as rj — F{<p)/N. Functional differentiation of 
F with respect to <fi shows that its jdobal maximum is at- 
tained for (f> — <p. Hence 77 ^ F((f>)/N for any choice of 
the test function (j). Choosing the exact solution (f> x of the 
phantom case yields a lower bound for the polymers with 
excluded volume 



1=1 



cq 



(11) 



where we used isotropy. 

Fixman [21] has also given an upper bound for the vis- 
cosity 

C N 

^jj^2(( K ^ R i-<) 2 ) ( 12 ) 



»=i 



cq 



in terms of N test vector fields {iij}, i = 1, . . . , N, subject 
to the divergence condition 



N 

E 



d 
dR~ 



{tpeqUi) = 0. 



(13) 



The problem with this estimate lies in the difficulty to find 
test fields that satisfy condition (13). The upper bound 
would match the lower bound (11) precisely, if we chose 
u\ — K V ^5R^ . However this choice in general does not 
fulfil condition (13) due to excluded-volume interactions 
between monomers in different clusters. The above test 



field is a valid choice for a system which consists just of 
a single cluster W. In that case the upper bound (12) 
coincides with the lower bound (11) and the equality 



.2 



(14) 



holds exactly. The additional subscript 's' of the thermal 
average indicates that there are no other degrees of free- 
dom besides the positions of the monomers in the cluster 
to be averaged — in contrast to the average in (11), which 
includes all monomers from all clusters. Eq. (14) has been 
suggested on a phenomcnological basis in the context of 
scaling arguments [8,11,25-31], which relate the viscosity 
to the longest time scale and the latter to the clusters' 
linear dimension i? gyr and the diffusion constant. 

Disorder average over crosslinks. — To obtain the 
viscosity of the macroscopic polymeric melt, we have to 
perform the macroscopic limit N — > 00 and average over 
different crosslink realisations. We focus on gelling sys- 
tems for which the appropriate crosslink ensemble is given 
by (i-dimensional lattice-bond percolation with bond prob- 
ability c/d, corresponding to the crosslink concentration c. 
We denote this combined procedure by an overbar • . For 
the viscosity of a polymeric melt, Eq. (11) yields the lower 
bound 

f 00 , s 

^Yd E^W^cq ( 15 ) 
n=l 

( n \ 

with f(n) denoting the cluster-size distribution and • 
the (normalised) partial average obtained from restricting 
• to clusters of fixed size n. Next we employ the scaling 
description of bond percolation for e := |c — c cr it| < 1 to 
compute critical exponents. In particular, the percolation 
average of the radius of gyration of n-clusters is known to 
obey the scaling relation [2] 



cq 



(n) 



,2/d f 



(16) 



with the fractal dimension df = {av)^ 1 w 2.53. It is ex- 
perimentally well confirmed (e.g. in [32,33]) and generally 
accepted [1,34] that the scaling (16) with the same fractal 
dimension also applies to gelation clusters with excluded- 
volume interactions in a melt. Moreover, on the theoreti- 
cal side, thermostatic properties of vulcanization clusters 
at criticality were shown to fall into the universality class 
of lattice-bond percolation [35]. Thus, when inserting (16) 
into (15), we obtain that the averaged viscosity of the poly- 
disperse melt diverges with a critical exponent 



k > 2v - (3 



(17) 



at the sol-gel transition. We recall that throughout our 
calculations we have not made any assumptions on the 
detailed nature of the excluded-volume potential. 

On the other hand, if excluded-volume interactions were 
absent, the geometric fractal dimension df would take on 
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(G) 

the Gaussian or ideal Rouse value d f « 4. The result- 
ing lower bound for k would then coincide with the ex- 
act result in the phantom case that was computed previ- 
ously [22-24]. 

Finally, we turn to the viscosity (14) of a single cluster 
with excluded-volume interactions that has been isolated 
from a gelling melt. On average, we get 



Vn 



c 



-(„) 



2d 



i eq, s 



(18) 



which is an exact result. For single clusters isolated from 
critical melts, (18) implies the Mark-Houwink relation 
Vn,s\e=o ~ n 2 / df s with the fractal dimension df tS for a sin- 
gle cluster. Whereas for clusters in the melt the excluded 
volume interaction is (partially) screened, this is not the 
case for an isolated cluster, resulting in a lower value of <if ;S 
as compared to df. A generalisation of Flory's argument 
to clusters [36] suggests df jS = 2. This exponent value is 
observed in light-scattering experiments for solutions ob- 
tained from diluting a polydisperse melt [16,37,38]. 

Discussion and outlook. — The exact lower bound 
for the critical exponent is definitely in disagreement with 
several entries in Table 1. For simulations or experiments 
referring to a melt, the only way out are effects due to en- 
tanglements (besides crossover effects due to a tiny critical 
region). Entanglement effects cannot be ruled out, but are 
generally not expected to play a vital role on the sol side 
of the transition quite in contrast to the well crosslinked 
gel phase, where entanglement effects due to interlocking 
loops are clearly important [39,40]. 

As far as dense solutions are concerned, one needs to 
consider hydrodynamic interactions due to the presence 
of a solvent. It seems feasible to include pre-averaged hy- 
drodynamic interactions in the variational approach pur- 
sued here. This may provide a semi-microscopic starting 
point for the computation of upper and lower bounds for 
viscosities for solutions of gelling polymers. Such bounds 
are all the more interesting, because a dense solution of 
gelling phantom polymers is known to behave completely 
unphysically, showing a finite viscosity at the sol-gel tran- 
sition [41]. 
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